arXiv: 1508.00002v2 [astro-ph.HE] 4 Aug 2015 


MNRAS 000, 1-10 (2015) 


Preprint 5 August 2015 


Compiled using MNRAS style file v3.0 


Constraining the parameters of the pulsar wind nebula 
DA 495 and its pulsar with Chandra and XMM-Newton 

A. Karpova^’^*, D. Zyuzin^, A. Danilenko^, Yu. Shibanov^’^ 

^loffe Institute, Politekhnicbeskaya 26, St. Petersburg, 194021, Russia 

^ Peter the Great St.Petersburg Polytechnic University, Politekhnicbeskaya 29, St. Petersburg, 195251, Russia 

Accepted XXX. Received YYY; in original form ZZZ 


ABSTRACT 

We present spectral and timing analyses of the X-ray emission from the pulsar wind 
nebula DA 495 and its central object, J1952.2-1-2925, suggested to be the pulsar, using 
archival Chandra and XMM-Newton data. J1952.2-1-2925 has a pure thermal spec¬ 
trum which is equally well fitted either by the blackbody model with a temperature 
of Ri 215 eV and an emitting area radius of Ri 0.6 km or by magnetized neutron star 
atmosphere models with temperatures of 80-90 eV. In the latter case the thermal 
emission can come from the entire neutron star surface which temperature is consis¬ 
tent with standard neutron star cooling scenarios. We place also an upper limit on the 
J1952.2-1-2925 nonthermal flux. The derived spectral parameters are generally com¬ 
patible with published ones based only on the Chandra data, but they are much more 
accurate due to the inclusion of XMM-Newton data. No pulsations were found and we 
placed an upper limit for the J1952.2-1-2925 pulsed emission fraction of 40 per cent. 
Utilizing the interstellar absorption-distance relation, we estimated the distance to 
DA 495, which can be as large as 5 kpc if J1952.2-1-2925 emission is described by 
the atmosphere models. We compiled possible multi-wavelength spectra of the neb¬ 
ula including radio data; they depend on the spectral model of the central object. 
Comparing the results with other pulsar plus wind nebula systems we set reasonable 
constraints on the J1952.2-1-2925 spin-down luminosity and age. We suggest that the 
Fermi source 3FGL J1951.6-1-2926 is the likely y-ray counterpart of J1952.2-1-2925. 

Key words: ISM: supernova remnants - ISM: individual: G65.7-I-1.2 - stars: neutron 
- stars: individual: J1952.2-1-2925 


1 INTRODUCTION 

More than 50 pulsar wind nebulae (PWNe) are currently 
known (e.g., Kargaltsev & Pavlov 2008). The Crab nebula 
is a classical example of young PWNe which is the best stud¬ 
ied from the radio to the TeV range. However, other PWNe, 
not yet studied in such details, show a variety of individual 
properties different from the Crab, which are determined by 
their pulsar parameters, evolution stages and surroundings, 
including the interstellar medium and the supernova ejecta. 
Studying such PWNe is crucial for understanding the prop¬ 
erties of their central ‘engines’, the interaction of the rela¬ 
tivistic pulsar wind with surroundings and particle accelera¬ 
tion mechanisms (Gaensler & Slane 2006). Multiwavelength 
observations of the PWNe at late evolutionary stages and 
different environment are still rare but they are particularly 
important for comparison with the Crab and various theo¬ 
retical models of evolution of these objects (Gelfand et al. 
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2009; Tanaka & Takahara 2010; Ghevalier & Reynolds 2011, 
and references therein). 

DA 495 (G65.7-I-1.2) appears to be an example 
of an evolved PWN. It was discovered in the radio 
(Galt & Kennedy 1968). Detailed radio observations showed 
that DA 495 has an uncommon annular morphology with 
a central flux deficit within -^2 arcmin from the centre, 
dubbed the radio ‘hole’ (Landecker & Gaswell 1983). The 
full extent of DA 495 is about 25 arcmin and the radio 
emissivity decreases gradually to its outer edge without any 
evidence of the supernova remnant (SNR) shell existence 
(see Fig. 1, bottom panel). A possible association with the 
open cluster NGG 6834 implied a distance D to the ob¬ 
ject of « 2.1 kpc (Landecker & Caswell 1983). Hi absorp¬ 
tion measurements and a flat rotation model for the Milky 
Way provided D within a range of 1-1.5 kpc, depending on 
the adopted Galactocentric radius for the Sun (Kothes et al. 
2004, 2008). Kothes et al. (2008) found a spectral break at 
1.3 GHz, the lowest for the known PWNe, with spectral in¬ 
dices a,, = 0.45T0.20 and a,y = 0.87T0.10 below and above 
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the break, respectively^. They also suggested that DA 495 
is about 20 kyr old. No detection of an associated pulsar in 
the radio and y-rays was reported. 

In X-rays, ROSAT and ASCA observations of DA 495 
revealed a faint compact source, IWGA J1952.2+2925 
(hereafter J1952), which is apparently surrounded by a dif¬ 
fuse non-thermal emission and is projected on the edge of the 
radio hole (Arzoumanian et al. 2004) (see Fig. 1). It was pro¬ 
posed to be a magnetospherically active neutron star (NS) 
powering the PWN, although only an upper limit on the 
pulsed emission fraction of 50 per cent for periods > 30 ms 
was derived from the analysis of the ASCA data. The com¬ 
pact source J1952 and its diffuse emission were later hrmly 
confirmed by Chandra high spatial resolution observations 
(Arzoumanian et al. 2008). It was established that J1952 is 
the point source located in the centre of the X-ray nebula 
with an extent of ~ 40". The nebula does not show any 
Crab-like torus+jet structure, however its spectrum is de¬ 
scribed by a power law with a photon index T = 1.6+0.3 typ¬ 
ical for PWNe. The latter allowed to state that it is the X-ray 
counterpart of the DA 495 PWN. J1952, presumably the pul¬ 
sar, has a pure thermal spectrum. It can be described either 
by the blackbody model with a temperature T « 2.5 MK 
and an emitting area radius R « 0.3 km or by the neutron 
star atmosphere (NSA) model (Zavlin, Pavlov & Shibanov 
1996) for a NS with the effective temperature of « 1 MK 
and the radius of 10 km. The interstellar absorption column 
density for the former and latter cases was « 2.3 x 10^^ and 
6.0 X 10^^ cm“^, respectively, which is a factor of 1.3-3.5 
lower than the entire Galactic absorption in this direction. 
Because of small count statistics, the spectral parameters 
were poorly constrained. In addition, two key parameters, 
the distance and the pulsar spin-down energy loss E, also re¬ 
mained uncertain, which did not allow to establish firmly the 
DA 495 PWN evolution stage. Finally, a high energy source 
3FGL J1951.6+2926 recently detected with Fermi/LAT was 
proposed as a possible y-ray counterpart (Acero et al. 2015) 
of the nebula. 

Here we report a simultaneous analysis of the Chan¬ 
dra^ and unpublished XMM-Newton^ X-ray archival data 
on DA 495. We also include in our analysis the extinction- 
distance relation towards J1952. This allows us to improve 
considerably the count statistics, to get an independent 
DA 495 distance estimate and to set more stringent con¬ 
straints on the PWN and pulsar parameters. We also use 
high temporal resolution XMM-Newton/EPlC-pn data to 
search for periodic pulsations from J1952 and derive a more 
stringent upper limit on its pulsed emission fraction. The 
details of observations are described in Sect. 2. The timing 
analysis is presented in Sect. 3. The extinction-distance re¬ 
lation and spectral analysis are described in Sect. 4 and 
5. In Sect. 6 we discuss our results, compare them with 
the DA 495 radio and y data and the data for other pul- 
sar+PWN systems. A summary is given in Sect. 7. 


^ Assuming the spectral flux oc where n is the fre¬ 

quency. 

^ PI Arzoumanian, Chandra/AClS-l, ObsID 3900 
® PI Arzoumanian, XMM-Newton/EPlC, ObsID 0406960101 
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Figure 1. 2.'5x2l5 X-ray Chandra {top), XMM-Newton {mid¬ 
dle) and 30'x27' radio VLA {bottom) images of DA 495 in 0.3-8 
keV range and at 1.42 GHz, respectively. The Chandra image was 
smoothed with a 3 pixel Gaussian kernel (logarithmic brightness 
scale is used). J1952 is clearly seen in the centre of the PWN in X- 
rays. The XMM-Newton image was binned to a pixel size of 1'.'6 
and smoothed with a 3 pixel Gaussian kernel (square root bright¬ 
ness scale is used). Solid and dashed circles depict the extraction 
region for the PWN+J1952 and background, respectively. The 
white circle in the radio image shows the position and extent of 
DA 495 in X-rays. 
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Figure 2. Data points show empirical Aiji-distance relation for 
the DA 495 direction derived using Marshall et al. (2006) extinc¬ 
tion map. The error bars represent Icr uncertainties. The solid 
line and gray filled region are smoothing spline approximations 
to the data points and their uncertainties, respectively. 


2 THE X-RAY DATA 

The XMM-Newton observations of DA 495 were carried on 
2007 April 21 with total exposure of about 50 ks. The EPIC- 
MOS cameras were operated in the Full Frame Mode with 
the medium filter setting, and the EPIC-pn camera was op¬ 
erated in the Small Window Mode with the thin filter. The 
XMM-SAS V. 13.5.0 software was used to process the data. We 
selected single and double pixel events (PATTERN < 4) 
for the EPIC-pn and single to quadruple-pixel events (PAT¬ 
TERN < 12) for the EPIC-MOS data. We removed periods 
of background flares using 10-12 keV and 12-14 keV light 
curves for the MOS and pn data, respectively. 

We also used the observation of DA 495 made with 
C/iandra/ACIS-I on 2002 December 9. The data mode was 
VFAINT, the exposure mode was TE and the total expo¬ 
sure was 24.8 ks. The ciAO v.4.6 chandra_repro tool with 
CALDB v.4.6.1 was used to reprocess the data. We found no 
background flares in this observation. The resulting effective 
exposures are presented in Table 1. 


3 TIMING ANALYSIS 

The EPIC-pn data, which have the time resolution of 5.6 
ms, were used to search for pulsations from J1952. Event 
times were corrected to the solar system barycentre using 
the SAS task barycen and J1952 Chandra coordinates ob¬ 
tained with the CIAO tool wavdetect (RA = 298?0710, Dec 
= 29?4313). We used a circular region of 20 arcsec around 
J1952 to extract its light curve. The SAS task epiclccorr 
was used for background subtraction and correction for vi¬ 
gnetting, bad pixels and quantum efficiency. The resulting 
light curve contained about 1060 source counts in the 0.5- 
1.4 keV energy band. We chose this energy range to mini¬ 
mize the contribution of the PWN emission and since J1952 
dominates at these energies (see Sect. 5). After that we per¬ 
formed a Fast Fourier Transform (FFT) analysis using stan¬ 


dard XRONOS v.5.22'‘ task powspec. No pulsations were 
found. We derived an upper limit for the pulsed fraction 
(PF) of 40 per cent (at the 99 per cent confidence level) 
for periods > 12 ms following the procedure described by 
Vaughan et al. (1994). 


4 DISTANCE AND INTERSTELLAR 
EXTINCTION 

It is possible to estimate the distance to DA 495 using 
the absorption column density An-distance relation along 
its line of sight. To do this we use the interstellar ex¬ 
tinction distribution obtained by Marshall et al. (2006) for 
the direction towards DA 495. They used the 2MASS sur¬ 
vey along with the BesanQon model of population synthesis 
(Robin et al. 2003) to calculate the Ks extinction depen¬ 
dence on distance along different lines of sight. We obtained 
this dependence from the VizieR Service® for the DA 495 
direction and then transformed the extinctions to the col¬ 
umn densities Ah utilizing a relation Ak^/Av « 0.1 (see 
e.g. Rieke & Lebofsky 1985) and a standard empirical Ah- 
Av relation (Predehl & Schmitt 1995). The resulting Ah- 
distance dependence is shown in Fig. 2. The An-distance 
relation and its uncertainties were approximated by smooth¬ 
ing splines. 


5 ANALYSIS OF THE X-RAY SPECTRA 

The Chandra and XMM-Newton images of DA 495 for the 
0.3-8 keV energy range are shown in the top and middle pan¬ 
els of Fig. 1, respectively. In the Chandra image the extended 
PWN is seen around the point source J1952. Due to a mod¬ 
erate angular resolution of XMM-Newton the PWN+J1952 
system appears to be slightly more extended® and the point 
source is blended with the PWN. Therefore, to perform spec¬ 
tral analysis of J1952 and the PWN we used their individual 
spectra extracted from the Chandra data, while from the 
XMM-Newton data we extracted the spectrum correspond¬ 
ing to the combined PWN-I-J1952 emission. For the latter we 
chose 30 arcsec aperture encircling the whole PWN emission 
(see Fig. 1) using evselect tool. We then used SAS tasks 
RMFGEN and ARFGEN to generate the redistribution matrix 
and ancillary response files for the MOSl, MOS2 and pn 
detectors. For the Chandra data set we extracted the J1952 
spectrum from a circular region of 1.5-pixel radius centred 
at this source and the PWN spectrum from the same cir¬ 
cular region as for XMM-Newton excluding 1 arcsec radius 
region around J1952 using SPECEXTRAGT tool. Background 
was taken from the region free from any sources (see Fig. 1). 
All spectra were grouped to ensure > 5 counts per energy 
bin. The resulting number of counts in the apertures and 
background-subtracted source count rates are presented in 
Table 1. 


http://heasarc.gsfc.nasa.gov/docs/xanadu/xronos/xronos.html 
® http://vizier.u-strasbg.fr/viz-bin/VizieR-3?- 
source=J/A%2bA/453/635/ 

® We have searched for a fainter PWN emission at larger spacial 
scales in the XMM-Newton images but not found it. 
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Table 1. Resulting Chandra and XMM-Newton exposure times, number of counts and background-subtracted count rates in 0.3-10 keV energy band for different instruments and 
source apertures for the J1952, PWN and J1952-I-PWN after data reprocessing, filtering and correction for background flares. 




Instrument 

Detector 

Chandra 

ACIS-I 

MOSl 

XMM-Newton 

EPIC/ 

MOS2 

pn 

Source 

J1952 

PWN 


J1952+PWN 


exposure (ks) 

24.8 

24.8 

46.5 

47.5 

33.7 

counts within the aperture (see text) 

145 

577 

1300 

1282 

4194 

background-subtracted count rate (cts ks“^) 

5.8±0.5 

19.6±1.0 

22.9±0.8 

22 . 2 ± 0.8 

74.4±2.2 


Table 2. Best-fit spectral parameters. All errors correspond to 90 per cent credible intervals. R and T are given as seen by a distant observer. The PWN luminosity Lpwn is in the 0.5—8 
keV range. A^pwn is the normalisation constant for the PL model of the PWN. -FJ 1952 is the point source bolometric flux. i® ^ 3cr upper limit on the non-thermal component of 

the point source luminosity in 0.5—8 keV band. 


Model 

A'h, 

D, 

R, 

T, 

P pwn 

-R^pwn 5 

LogFji952, 

LOgLj,g 52 , 

LogZ/pwn; 

X^/dof 


10^^ cm“^ 

kpc 

km 

eV 


10-5 ph 

keV“^ cm“3 

erg s“^ cm“^ 

erg s“^ 

erg s“^ 


1. BB 

2 e-l-0.5 


'-'•^-0.3 

215^23 

1 yi +0.12 
-*-• ‘ -*--0.12 

4 4+0-'^ 

^•^-0.6 

iz.y^_0 11 

< 31.51 

q9 99+0.38 

455/461 

2. NSMAX 

B = IOI 2 G 

N 0 

d d 

+ 1 

CO 

0 0 + 1-7 

'^•'^-1.3 

lot^i 


1 oq+0.14 
-*--°'^-0.13 

c 1 +0.9 
^•-‘■-0.7 

— 12 47+'^-35 

< 31.66 

q9 cq+0.37 

OZ.OO_Q 45 

455/461 

3. NSMAX 

B = I 013 G 

0 4+0.7 

q q + 1-5 
'^•^-1.2 

(^+9 

’^-3.5 

9AI5 

1 09 + 0.13 
t-o^_0,l3 

^ 1+0.8 
^•-*--0.6 

— 12 ^^7+0.27 

iz.o /_0 20 

< 31.72 

32.50l°i3 

454/461 


to 

o 

or 
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A power law (PL) was used to model the PWN emis¬ 
sion. For the J1952 spectrum, we tried the neutron star 
atmosphere models NSA (Zavlin et al. 1996) and NSMAX 
(Ho, Potekhin & Chabrier 2008) with a gravitational red- 
shift Zg + 1 = 1.21 and magnetic fields B = 10^^ G and 
B = 10^® G and the blackbody (BB) model. For the in¬ 
terstellar absorption, we used the xspec photoelectric ab¬ 
sorption PHABS model with default cross-sections bcmc 
(Balucinska-Ghurch & McCammon 1992) and abundances 
ANGR (Anders & Grevesse 1989). We fitted XMM-Newton 
and Chandra spectra of the J1952+PWN, J1952 and PWN 
simultaneously in 0.3-10 keV energy range. 

We fit the data using the Bayesian inference with 
the Nn-distance relation as a prior. We obtained samples 
from the posterior distribution by the Markov chain Monte- 
Garlo (mcmc) sampler developed by Goodman & Weare 
( 2010 ) and implemented as a python package emcee by 
Foreman-Mackey et al. (2013). We obtained 100000 sam¬ 
ples which is enough compared with autocorrelation time 
(50-80). The procedure is similar to that utilised by 
Kirichenko et al. (2015). Posterior median values of the 
spectral parameters with 90 per cent credible intervals are 
presented in Table 2, where the models differ only by the 
thermal components applied and are named according to 
that. Since the best-fit parameters do not depend strongly 
on the specific type of atmospheric models, i.e., NSA or NS¬ 
MAX, we present the results only for the BB and NSMAX 
models. The values per degree of freedom (dof) show 
that all models describe the data equally well. For both 
models our results are generally consistent with results of 
Arzoumanian et al. (2008), however they are several times 
less uncertain (cf. their Table 1) due to the better count 
statistics. The example of the fit for the model 2 from Table 
2 is shown in Fig. 3. 

It is important, that using the Bayesian approach and 
Nb_-D relation, we were able to derive the values of D and 
the J1952 emitting area radius R independently, as opposed 
to their ratios obtained traditionally from such fits. For the 
NSMAX models, R = IQff km (B = lO^^ G) and R = 
big g km [B = 10^® G) imply that the thermal emission 
originates from the entire NS surface. For the BB model, 
the emitting area is significantly smaller, R = O.OIq ® km, 
and the radiation likely comes from a hot polar cap of the 
pulsar. To derive an upper limit on the temperature of the 
entire NS surface in the latter case, we added a second BB 
component to the best-fit model. The radius of the emitting 
area was fixed at a conventional NS radius of 13 km. We 
obtained a 3a upper limit on the entire surface temperature 
of « 60 eV. 

For the J1952 spectrum we also tried the PL model that 
provided a worse fit quality (x^/dof=489/461) and resulted 
in a too large photon index F = not typical for 

pulsars (see e.g. Li, Lu & Li 2008). 

Besides, we added a power law component to the ther¬ 
mal component for J1952 to estimate an upper limit on the 
pulsar non-thermal luminosity Ljlg 52 - In current data this 
emission component is not resolved within the non-thermal 
PWN background. The photon index F was assumed to lie 
in the range 0.5 < F < 2.0 typical for other pulsars. The 
resulting values of ^^1952 are also presented in Table 2. 


6 DISCUSSION 

We used high temporal resolution XMM-Newton/EPlC-pn 
archival data to search for pulsations from the NS J1952. 
No pulsations were detected. However, we set a more strin¬ 
gent upper limit for the pulsed fraction of 40 per cent 
(99 per cent confidence) in a wider period range of > 12 ms 
(Sect. 3) than previous estimates based on the ASCA 
data (Arzoumanian et al. 2004). We confirmed the result of 
Arzoumanian et al. (2008) that only the thermal emission 
component is detected in the X-ray spectrum of J1952 and 
better constrained its parameters. For the BB model the 
emitting area radius is « 0.6 km (see Table 2) showing that 
the thermal component originates from a small hot spot on 
the NS surface with T « 215 eV, which is typical for pul¬ 
sar hot polar caps. On the other hand, the NSMAX models 
suggest that the emission comes from the entire NS surface. 
We cannot assert definitely which case is true for J1952. 
The detection of the pulsations, the phase-resolved spectral 
analysis and the analysis of pulse profile shapes could help to 
solve this problem. Our upper limit on the PF of 40 per cent 
is higher than typical PFs of 10-30 per cent detected for the 
thermal emission from entire surfaces of other NSs (Zavlin 
2007; Pavlov & Zavlin 2000). 

The interstellar absorption-distance relation (Sect. 4) 
allowed us to estimate the distance to DA 495. For the 
BB model D = 2Af_\ l kpc. The NS atmosphere models 
suggest a larger distance: kpc, if the J1952 mag¬ 

netic field is 10^^ G, or 3.2'f\'2 kpc, if it is lO*^® G (Ta¬ 
ble 2). Within uncertainties, all these values are consistent 
with one of the previous estimates, D « 2.1 kpc, assum¬ 
ing possible association of DA 495 with the open cluster 
NGC 6834 (Landecker & Gaswell 1983). Another distance 
estimate of 1.0 ± 0.4 kpc is compatible only with the BB 
model. It was based on the Hi absorption feature in the 
DA 495 polarized radio continuum emission with a systemic 
velocity of - 1-12 km s“^, which is likely associated with the 
Galactic Local Arm (Kothes et al. 2008). The authors also 
mentioned a distance of 5 kpc as an alternative interpreta¬ 
tion of the feature. This is compatible at the upper confi¬ 
dence level with the result provided by our atmospheric fits. 
However, they argued that this distance is inconsistent with 
the low foreground absorbing column density obtained by 
Arzoumanian et al. (2008) from Chandra X-ray data using 
the BB model and with the absence of the Hi absorption 
from the Local Arm tangent point expected to be at the 
velocity of about +20 km s“^. We note, that the latter ab¬ 
sorption cannot be completely excluded due to a low signal- 
to-noise ratio of the radio data. The atmospheric models 
cannot be also ignored in the above argument. Moreover, re¬ 
cent trigonometric parallax measurements of high mass star 
forming regions show that the Local Arm extends beyond 
5 kpc with a shallow pitch angle of « 10° between I « 72° 
and 52° toward the Perseus Arm (Reid et al. 2014; Xu et al. 
2013). It is possibly a separate branch of the Perseus, con¬ 
necting with it at a distance of about 6.5 kpc (see fig. 12 
from Xu et al. 2013). DA 495 with I « 65? 7 is obviously lo¬ 
cated within the new Local Arm longitude range. Therefore, 
the gas within this arm can be a natural source of the larger 
foreground X-ray absorption at the larger DA 495 distance 
suggested by the atmospheric models (Table 2). 

The PWN spectral indices derived from our X-ray 
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Figure 3. Observed X-ray spectrum of the DA 495 pulsar-|-PWN system obtained with XMM-Newton (top-left) and Chandra (top- 
right). The data for different instruments/sources are shown by different crosses, as indicated in the inserts. Lines show the best-fit 
model. Bottom: Fit residuals. 



Figure 4. Unabsorbed spectrum of the DA 495 PWN from the 
radio to X-rays. The X-ray spectrum of the associated pulsar 
J1957 is fitted by the BB model. The dashed and solid lines are 
the best fits to the PWN data by the broken PL model with the 
single and the double-knee breaks, respectively. The position of 
the second high frequency spectral break is poorly defined (see 
text for details), it is fixed here for illustrative purposes. 

fits are 0.71 ± 0.12, 0.83 ± 0.13 and 0.82 ± 0.13 (90 per cent 
confidence) when the thermal component from the NS was 
described by the BB and atmosphere models with B = 
10^^ G and 10^® G, respectively. These values are in agree¬ 
ment within uncertainties with each other and with a,, = 
0.87 ±0.10 (68 per cent confidence) obtained from the radio 
data for frequencies above the 1.3 GHz break (Kothes et al. 
2008). However, a marginal spectral flattening in X-rays, 
particularly in the BB case, cannot be excluded, suggesting 
a second break. To examine whether an additional break be¬ 
tween the radio and X-rays is needed or not, we tried to fit 
the X-ray spectra together with the radio data taken from 
Kothes et al. (2008). At first, the PWN spectrum was fitted 
by a broken power law (BPL) with a single break. In all 


cases the break in the radio was 2.lj)Q g GHz. This value is 
higher than that of Kothes et al. (2008) due to the addition 
of the X-ray data; using only the radio data we reproduced 
their result of « 1.3 GHz. The spectral indices and q „2 
before and above this break were 0.61 ±0.07 and 1.00 ±0.01 
for the NSMAX±BPL models and 0.58±0.07 and 1.00±0.01 
for the BB±BPL model. Then we included a second break 
to the PWN model. For the NSMAX±BPL models, the re¬ 
sulting spectral indices Oj/i and q „2 before and above the 
break at 2 GHz and the spectral index a „3 above the sec¬ 
ond break were O.SOIq ° 2 , 1 - 021 q q 4 and 0.84 ))q )2 in the case 
of B = 10^^ G and 0.50l^;)^, 1.02t°;J° and 0 . 831 ^ 1 ^ in the 
case of B = 10^^ G. For the BB±BPL model, a^i — 0.5^q °® , 
ctu 2 = I-OSIq (,5 and a „3 = 0.72tQ 44 . We carried out an F- 
test using the xspec ftest routine to compare models with 
single and double-knee breaks. The resulted probabilities for 
the NSMAX±BPL were « 0.1 for both values of magnetic 
field demonstrating that the second break is not needed. On 
the other hand, the probability for the BB±BPL model is 
0.03 providing a marginal evidence that the PWN spectrum 
in X-rays indeed becomes flatter although the second break 
position is poorly defined and lies in a range of 10 ^"* — 10 ^^ 
Hz, somewhere between the mid-infrared and X-rays. The 
best fits to the data by the BB±BPL model with single and 
double-knee breaks are presented in Fig. 4. 

Assuming that the spectral break at « 2 GHz in DA 495 
emission arises from the synchrotron cooling and using the 
formula from Ginzburg & Syrovatskii (1965), 

^c(GHz) « 1.2 X IO^B-^gC, (1) 

where i/c is the cooling frequency, BmG is the PWN mag¬ 
netic field in milligauss, and tkyr is the PWN age in kyr, 
we can estimate the age of DA 495. The synchrotron life¬ 
time is Tsynch « 3.3Bpcc/u yr, where Rpc is the emitting 
area radius in pc, v is the velocity of the pulsar wind parti¬ 
cles and c is the speed of light. At a distance of about 1.5 
kpc, which is close to the low boundary obtained from the 
X-ray spectral fits, Rpc = O.lfiBi.skpc, which corresponds 
to the 20 arcsec angular radius of the X-ray PWN, and 
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''"synch ~ 1-7-Di.skpc yr, assuming a typical velocity v = 0.3c 
(Kennel & Coroniti 1984). On the other hand, at an X-ray 
photon energy E = 1.5 keV, r^ynch = 1.2 £^kev^-®mG ^ ~ 
This yields BmO ~ 0-7£l)~5kpc expression (1) 

implies the DA 495 age tkyr ~ dODi.skpc. At the maximum 
distance of about 5 kpc we obtain SmG ~ 0.3, which is com¬ 
parable to that of the Crab Nebula, and a maximum age 
tkyr ~ 130. These simplified estimates imply that DA 495 
may be a factor of 2-6 older than it was suggested by 
Kothes et al. (2008). 

The distance estimate allowed us to derive the PWN 
luminosity Lpwn and the 3a upper limit on the non- 
thermal luminosity of its pulsar Tjlg 52 (Table 2). These 
parameters can be used to estimate the J1952 spin- 
down luminosity E utilizing the empirical dependencies 
Tpwn(f5), L^l^{E) and Lp^iL^ld obtained for 0.5-8 keV 
band by Kargaltsev & Pavlov (2008) using the data for 
other pulsars-fPWN systems (Fig. 5). For completeness, 
we added to these data PSR J2022-I-3842 with its PWN 
G76.9-1-1.0 (Arzoumanian et al. 2011; Arumugasamy et al. 
2014) marked in Fig. 5 by a triangle. Summarising the re¬ 
sults for all three models presented in Table 2, logI/p„n for 
DA 485 lies in the range of 31.7-32.9 (horizontal solid lines in 
the top and bottom panels of Fig. 5), which implies logE = 
32.7-37.9 for the X-ray efficiency range = 10“® — 10“^ 
based on the data for all PWNe observed (the grey shaded 
region between the horizontal solid lines in the top panel). 
Using the respective empirical luminosity vs age dependence 
Tpwn(t) for the same PWNe set (Fig. 6, left panel)^ and 
the DA 495 Lpwn allowed range, we can obtain the pos¬ 
sible age range of DA 495 as 7<fkyr^l55 (the grey re¬ 
gion between two vertical dashed lines in the left panel of 
Fig. 6)®. This age range is slightly wider than that based on 
the simple synchrotron cooling scenario. At the same time, 
the derived upper limit on the pulsar nonthermal luminos¬ 
ity logLjJggg < 31.7 (tho hortzontul solid line in the middle 
panel of Fig. 5) suggests logiJ < 37.1 for the pulsar efficiency 
range > 10“®'^ based on the data for all pulsars (the grey 
region in the middle panel of Fig. 5). The lower limit for ri^'^ 
(the dotted line in the middle panel of Fig. 5) corresponds to 
the Vela pulsar X-ray efficiency, which is known to be very 
under-luminous in X-rays. Consequently, a real logLl value 
likely lies in the range of 32.7-37.1 where the allowed ranges 
of E for the PWN and J1952 overlap. This range boundaries 
are shown by solid lines in the E{t) empirical dependence in 
the right panel of Fig. 6. On the other hand, the Lpwn(Lpsr) 
dependence allows us to set a lower limit of « 30.7 on the 
J1952 non-thermal luminosity (the vertical dashed line in 
the bottom panel of Fig. 5). This value is also shown by the 
horizontal dashed line in the middle panel of Fig. 5. Finally, 
the E{t) empirical dependence in the right panel of Fig. 6 
and the estimated age range (the vertical dashed lines in 


^ For age t, the pulsar spin-down ages are used except for 
the cases when the true ages are known from respective 
SNR studies: Crab, J1833—1034 (Bietenholz & Bartel 2008), 
J1811—1925 (Roberts et al. 2003), B1951-I-32 (Migliazzo et al. 
2002), J0538-I-2817 (Ng et al. 2007). 

® At this point we excluded PSR J1119—6127, marked by the 
star symbol in the left panel of Fig. 6, as a peculiar high-B pulsar 
with a very faint PWN (Gonzalez & Safi-Harb 2003). 
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Figure 5. Top panel: X-ray luminosities of PWNe vs the spin- 
down luminosity E. Middle panel: the same as is in the top panel 
but for X-ray non-thermal luminosities of the pulsars powering 
the PWNe. Bottom panel: X-ray luminosities of the same PWNe 
vs pulsar X-ray non-thermal luminosity. Error-bars include 40 per 
cent systematic distance uncertainties for objects where parallax 
was not measured. Grey solid lines show constant efficiency lev¬ 
els r) for the PWNe and PSRs and their ratio as indicated in the 
plots. Horizontal solid lines correspond to the DA 495 PWN lumi¬ 
nosity lower and upper bounds and the upper limit on the J1952 
pulsar non-thermal luminosity obtained from the X-ray spectral 
fits. Their values are presented in the plots. The dotted line in 
the middle panel indicates the lowest efficiency level correspond¬ 
ing to the Vela pulsar. The black dashed lines represents the lower 
bound of the J1952 pulsar non-thermal luminosity. The regions 
where the DA 495 PWN and its pulsar are expected to locate are 
filled by the grey colour. 
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Figure 6 . Lpwn for the same PWNe as in Fig. 5 (left) and respective E (right) vs age t. Solid lines show the DA 495 PWN luminosity 
bounds derived from the fits and its spin-down luminosity range estimated from Lpwn(F) and Lp]„.(E) dependencies presented in Fig. 5. 
The vertical dashed lines show a possible age bounds for DA 495 and the horizontal dashed line in the right panel demonstrates its 
possible spin-down luminosity lower bound based on the age bounds and E(t) distribution for other pulsars. The regions where the 
DA 495 PWN and its pulsar are expected to locate are filled by the grey colour. 


this plot) suggest a more narrow \ogE range of 34.7-37.1 
corresponding to the grey region in this panel between the 
upper solid and horizontal dashed lines. This appears to be 
the most reasonable DA 495 E range summarising all data 
available for DA 495 and other PWN-|-pulsar systems. 

Comparing thermal properties of J1952 derived from 
the X-ray spectral fits (Table 2) with NS cooling models 
(Yakovlev & Pethick 2004), we find that the entire surface 
temperature of 80-90 eV provided by the atmospheric mod¬ 
els is consistent with the standard cooling scenario for the 
obtained age range. For the BB model, a 3cr upper limit 
on the entire surface temperature of «60 eV appears to be 
cooler than the standard cooling scenario predicts, if J1952 
age is < 30 kyr. 

Arzoumanian et al. (2011) have noticed, that by many 
properties DA 495 is similar to the PWN G76.9-1-1.0 associ¬ 
ated with the pulsar J2022-I-3842. Both objects demonstrate 
bright radio PWNe with similar morphologies, while their X- 
ray PWNe are faint, have much smaller sizes than in the ra¬ 
dio and coincide with central faint emission holes of the radio 
PWNe. The ratio of the PWN sizes in the radio and X-rays is 
about 20 for both objects. As for DA 495, no SNR shell was 
found for G76.9-bl.O. PSR J2022+3842 with P = 48.6 ms, 
detected in the radio and X-rays, has a characteristic age 
« 8.9 kyr, B « 2.1 x 10^^ G and B « 3 x 10®'^ erg s"^ 
Its PWN X-ray luminosity is « 6 x 10^^ erg s“^, as¬ 
suming the distance of 10 kpc (Arzoumanian et al. 2011; 
Arumugasamy et al. 2014). Positions of PSR J2022-I-3842 
and its PWN in Fig. 5, 6 are shown by triangles. Given the 
similarity between DA 495 and G76.9-1-1.0, one can suggest 
that DA 495 may also be as young as G76.9-1-1.0, which is 
consistent with the low boundary of the age range predicted 
by the E — Lpwn relation (see Fig. 6 ). For the simple syn¬ 
chrotron cooling model the minimum age is about 30 kyr at 
the minimum allowed distance of 1.3 kpc inferred from the 
X-ray fits. This demonstrates that the simple synchrotron 
cooling model does not work in this case. This is not a sur¬ 


prise, as it cannot explain spectra of other PWNe also and is 
used only for rough estimates (e.g. Gaensler & Slane 2006). 

Fermi/LAT recently detected at 6 cr significance a 7 -ray 
source 3FGL J1951.6-1-2926 coinciding spatially with DA 495 
(Acero et al. 2015). Its spectrum is described by a power law 
with an exponential cut-off with an integral flux « 10 “^^ 
erg cm“^ s“^ above 100 MeV. It was proposed as a possi¬ 
ble DA 495 PWN 7 -ray counterpart. However, comparing 
the source spectrum with the radio-X-ray spectrum of DA 
495 compiled by us above, we find that such interpretation is 
hardly possible, since it would yield an unrealistically strong 
high energy excess in the PWN spectral energy distribu¬ 
tion (SED) (Fig. 7) which is not predicted by any model of 
the PWN spectral evolution and not observed for most of 
evolved PWNe (e.g., Temim et al. 2013, 2015). 

An exception is a composite SNR MSH 11-62 contain¬ 
ing the PWN where similar excess is observed (Slane et al. 
2012). To explain the excess with SNR models one needs 
to assume an unrealistically high ambient density for a 
hadronic scenario or an unreasonably large fraction of the 
SNR mechanical energy converted to relativistic particles 
for a leptonic scenario. The PWN model, in turn, requires 
an extremely high pulsar spin-down power. Based on that, 
Slane et al. (2012) suggested that such a strong 7 -ray emis¬ 
sion of the SNR likely arises from its pulsar. 

Another unusual example is G76.9-1-1.0, the twin of 
DA 495 discussed above. Martin et al. (2014) modelled its 
spectrum and concluded that the detection of this PWN 
at the GeV-TeV range seems unexpected. Nevertheless, 
Acero et al. (2015) have reported on a possible 7 -ray coun¬ 
terpart, 3FGL J2022.2-b3840, of the G76.9-f 1.0 PWN. We 
find, that this implies similar unreasonably high energy ex¬ 
cess in its SED as for DA 495 and MSH 11-62 (Fig. 8 ). 
At the same time, Pilia et al. (2011) have reported on a 5a 
detection of the pulsed emission from PSR J2022-I-3842, as¬ 
sociated with this PWN, at energies above 200 MeV with a 
pulsed flux of (2.7±0.5) X10“^ photons cm“^ s“^. This flux 
is compatible with that reported by Acero et al. (2015) for 
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Figure 7. Radio and X-ray SED of the DA 495 PWN 
and the GeV spectrum of the associated Fermi source 3FGL 
J1951.6+2926. 



Figure 8. Radio and X-ray SED of the G76.9+1.0 PWN and the 
GeV spectrum of 3FGL J2022.2+3840, presumably the associated 
pulsar J2022+3842. The X-ray spectrum of the pulsar, marked by 
PSR, is shown for comparison. The y-ray emission level in respect 
to the radio-X-ray SED is similar to that of DA 495. 


3FGL J2022.2+3840, suggesting that both groups observed 
the same object. For the distance of 10 kpc, this yields the 
pulsar luminosity L.y = (3.6 ± 0.6) x 10^® erg s“^ and the 
y-efiiciency L-y/E « 0.01. 3FGL J2022.2+3840 has also the 
power law spectrum with the exponential cut-off typical for 
other Fermi pulsars (Abdo et al. 2013). 

The DA 495 SED similarity with those of G76.9+1.0 
and MSH 11-62 implies that the associated Fermi source 
3FGL J1951.6+2926 is likely the y-ray counterpart of the 
J1952 pulsar but not of the PWN. This is also supported 
by the pulsar-like spectrum of the source. In this case, the 
pulsar luminosity Ly « 10®^ erg s“^ for the distance of 3 kpc. 
Given that, and assuming the same y-efiiciency as for PSR 
J2022+3842, we obtain E ~ 10^® erg s“^ for J1952, which 
is within the range obtained above using the X-ray data. 


7 SUMMARY 

The X-ray spectrum of J1952 is pure thermal and can be de¬ 
scribed either by the BB or NS atmosphere models. In the 
former case the emission comes from a pulsar hot spot and 
in the latter - from a bulk of the NS surface with T «85 eV 
consistent with standard NS cooling scenarios. The Nn^D 
relation provides independent constraints on the distance to 
the object. For atmosphere models it can be larger, up to 
5 kpc, than it was supposed early. The upper limit for the 
pulsed emission fraction of J1952 is 40 per cent at 99 per cent 
confidence. It does not help to discriminate whether the 
thermal emission originates from the hot spot or the bulk of 
the NS surface. The derived PWN X-ray luminosity range 
and the pulsar nonthermal luminosity upper limit combined 
with those of other pulsar+PWN systems imply the DA 495 
age of 7-155 kyr. Empirical relations between E, Lpwn, Tpir 
and age predict a possible J1952 logE range of 34.7-37.1. 
The BB model for J1952 suggests that the PWN spectrum 
likely becomes flatter in X-rays than in the radio above a 
break near 2 GHz, suggesting a second high frequency spec¬ 
tral break. The atmosphere models do not require the second 
break. Sub-millimetre and infrared observations could help 
to check that. We suggest, that 3FGL J1951.6+2926 is the 
likely y-ray counterpart of the J1952 pulsar. To constrain 
J1952 energetics and age, deeper X- and y-ray timing obser¬ 
vations are needed. 
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